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ABSTRACT 

Context. High-energy observations have unveiled peculiar classes of isolated neutron stars which, at variance with radio pulsars, 
are mostly radio silent and not powered by the star rotation. Among these objects are the magnetars, hyper-magnetized neutron 
stars characterized by transient X-ray/7-ray emission, and neutron stars with purely thermal, and in most cases stationary, 
X-ray emission (a.k.a., X-ray dim isolated neutron stars or XDINSs). While apparently dissimilar in their high-energy behavior 
and age, both magnetars and XDINSs have similar periods and unusually high magnetic fields. This suggests a tantalizing 
scenario where the former evolve into the latter. 

Aims. Discovering so far uninvestigated similarities between the multi-wavelength properties of these two classes would be a 
further step forward to establish an evolutionary scenario. A most promising channels is the near infrared (NIR) one, where 
magnetars are characterized by a distinctive spectral flattening with respect to the extrapolation of the soft X-ray spectrum. 
Methods. We observed the two XDINSs RX J0420.0-5022 and RX J1856. 5-3754 with the Multi-Conjugate Adaptive Optics 
Demonstrator (MAD) at the Very Large Telescope (VLT), as part of the instrument guaranteed time observations program, to 
search for their NIR counterparts. 

Results. Both RX J1856.5-3754 and RX J0420.0-5022 were not detected down to K s ~ 20 and K s ~ 21.5, respectively. 
Conclusions. In order to constrain the relation between XDINSs and magnetars it would be of importance to perform deeper 
NIR observations. A good candidate is 1RXS J214303. 7+065419 which is the XDINS with the highest inferred magnetic field. 

Key words. Infrared: stars; neutron stars: individual RX J0420.0-5022,RX J1856. 5-3754 



1. Introduction 

Over the last two decades, high-energy observations have 
unveiled the existence of peculiar classes of isolated neu- 
tron stars (INSs) like the anomalous X-ray pulsars and 
the soft 7- repeaters (AXPs and SGR s), i.e., the magne- 
tar candidates (e.g., 
INSs (XDINSs; e.g., 



Mereghettill2008l) . 

Haber]||2007l) . Despite their vastly dif- 
ferent observational manifestations, both magnetars and 
XDINSs are slow rotators (P ~ 1-10 s), with surface mag- 
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netic fields reaching ss 10 14 -10 15 G in the magnetars case, 
have persistent X-ray luminosities much larger than the 
inferred rotational energy losses Lx ~ 10-100 E, and they 
are usually radio-silent. 

The XDINSs X-ray spectra are purely thermal and 
well represented by a blackbody (kT w 50-100 eV) whose 
emission radius is consistent with a sizable fraction of 
the neutron star surface. In most cases, one or possibly 
more broad absorption features (E# ne « 0.2-0.7 keV) have 
been observed (see Haberl 2007 and references therein), 
which hint toward magnetic fields of ~ 10 13 — 10 14 G, 
comparable to those derived from the spin down. The 
absence of power-law tails is consistent with the faint 
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rotation-powered magnetospheric emission expected from 
the inferred spin-down luminosity (E « 3 — 5 x 10 30 
erg s _1 ). In the optical, only a few XDINSs have been 
identified (see Kaplan 2008 and references therein) and 
a candidate counterpart was recently proposed for 1RXS 
J214303. 7+065419 (Zane et al. 2008), while that of RX 
J0420. 0—5022 was not confirmed (Mignani et al., in prepa- 
ration). Their spectra mostly follow a Rayleigh- Jeans, 
with fluxes exceeding by a factor > 10 the low-energy 
extrapo l ation of the X-rays best-fitting blackbody (e.g., 
Kaplan . 2008). Whcthci this is indeed due to emission 
from regions of the sta r surface at different temperature 
(e.g.. iPons et all 120021 ) is still under debate. Like in the 
X-rays, no rotation-powered magnetospheric emission is 
detectable in the optical, unless XDINSs are a factor of 
> 10 3 more efficient emitters than rotation-powered neu- 
tron stars (Zharikov et al. 2006). 

Whether there is a (evolutionary?) link between mag- 
netars and XDINSs is still hotly debated (see McLaughlin 
et al. 2006; Popov et al. 2006, for a discussion). As pro- 
posed by Mignani et al. (2007a), the detection of the op- 
tical/near infrared (NIR) flattening observed in the mag- 
netar spectra, ascribed either to magnetospheric emission 
powered by the magnetic field or to thermal emission from 
a fallback disk, would strengthen such a link. For the 
XDINSs, such flattening would be easily recognizable in 
the NIR, where the contribution of the optical Rayleigh- 
Jeans continuum is negligible. So far, NIR observations of 
XDINSs performed with the Very Large Telescope ( VLT) 
did not unveil counterparts down to H= 21.5-22.9 (Lo 
Curto et al. 2007). 

We report on new NIR observations of RX 
J1856. 5-3754 and RX J0420.0-5022 recently performed 
in the K s band with the VLT, the first presented for these 
sources. Observations are described in Sect. 2, while re- 
sults are presented and discussed in Sects. 3 and 4, re- 
spectively. 



2. Observations 

2.1. Observations description 

We performed K s -band observations of RX J1856.5— 3754 
and RX J0420.0-5022 on September 28 and 29 
2007 (MJD=54371 and 54372), respectively, using the 
European Southern Observatory (ESO) Multi- Conjugate 
Adaptive Optics Demonstrator (MAD)}], mounted at the 
Very Large Telescope ( VLT) Melipal telescope, as part 
of the guaranteed time program. The MAD instrument 
was developed within the European Extremely Large 
Telescope framework. The instrument accommodates two 
wavefront sensors (WFS): a star oriented (SO) multi 
Shack-Hartmann and a layer oriented (LO) multi pyra- 
mid (Regazzoni 2000; Regazzoni et al. 1996; 2000). Both 
sensors can use reference stars identified in the instrument 
2' field of view. The MAD instrument is completed by the 
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Camera for Multiconjugate Adaptive Optics (CAMCAO), 
a NIR camera equipped with a 2K x 2K Hawaii II infrared 
detector movable over the 2' corrected circular field-of- 
view. The pixel-scale is 0'.'028, corresponding to a 57" x 57" 
field-of-view on the detector. Thanks to its flexibility, the 
instrument can be used in classical single reference adap- 
tive optics (AO), ground layer adaptive optics (GLAO), 
or multi-conjugate adaptive optics (MCAO) according to 
the availability of suitable reference stars close to pointing 
direction. 

During the first night, we used the LO WFS in clas- 
sical AO mode, with a V w 11.3 reference star. During 
the second night, we used the MCAO mode with 4 fainter 
stars as a reference (V = 14.9, 15.3, 15.4 and 15.6). The 
seeing conditions during the observations were 1"23 for 
RX J0420. 0-5022 and l'/37 for RX J1856. 5-3754 and 
the airmass was 1.11 and 1.35, respectively. For RX 
J0420. 0—5022, we obtained a total of dithered 65 expo- 
sures along a pattern that covered a region of ~ 25" x 25". 
The exposures were taken in blocks of 5, i.e., along each 
node of the dithering pattern, with 30 s detector inte- 
gration time (DIT) and NDIT=2. The total integration 
time was 3900 s. WE performed the observation of RX 
J1856.5— 3754 in a similar way, with a dithering pattern 
covering a region of ~ 35" x 15". A total of 25 exposures 
were obtained with DIT=30 s and NDIT=4, for a total 
integration time of 3000 s. 

2.2. Data reduction and calibration 

For each set of science images we applied standard data 
reduction procedures for NIR observations, which include 
trimming of the frames, bad pixel masking, dark and flat- 
field correction, and sky-background subtraction. Special 
care was taken for images alignment, allowing a correc- 
tion for rotation to account for a non-optimal efficiency of 
the de-rotator. For the flat-field correction, we used the 
median of several images obtained on the sky at the be- 
ginning of each night. The bad pixel mask was created by 
analyzing the ratio of two flat -field images with signifi- 
cantly different exposure levels. For each science observa- 
tion we then constructed a reference sky image. Then we 
constructed the reference sky image from the median of 
all science images of a given observation. In the case of our 
targets, the fields are not affected by crowding of stellar 
sources, which allows for the creation of a good quality sky 
image. This was then normalized to the median counts of 
each science image and subtracted. Finally, we combined 
all sky-subtracted science images after properly aligning 
them using as a reference the positions of the brightest 
(not saturated) sources detected in all images. 

From the final co-added images we measured the as- 
sociated image quality from the average full width at half 
maximum (FWHM) of the intensity profile of a number 
of reference stars. Although the nominal seeing condi- 
tions in the night were not optimal (see above), thanks 
to the MAD AO correction, the FWHM of stellar im- 
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Fig. 1. We present 30" x 30" MAD K s -band images of the RX J1856. 5-3754 (left) and RX J0420.0-5022 (right) 
fields. North to the top, East to the left. The circles mark the expected target positions, with the radii (0"35 and 1", 
respectively) corresponding to their positional uncertainty (see Sect. 3.1). The labeling in the RX J1856.5— 3754 field 
is the same used in van Kerkwijk & Kulkarni (2001). The few faint objects detected in the RX J0420. 0—5022 field are 
also labeled. The intensity scale in the image has been adjusted for a better visualization. The edge of the bright field 
galaxy ESO 202-8 is visible in the top right of the image. 



ages measured on the detector actually turned out to be 
0"28 and 0"40 for the RX J0420.0-5022 and for the RX 
J1856.5— 3754 field, respectively. We performed the photo- 
metric calibration by observing NIR standard stars during 
the same nights. We cross-checked our calibration by di- 
rect comparison with objects of the 2MASS point-source 
catalog (Skrutskie et al. 2006). 

3. Data analysis and results 

3.1. Astrometry 

The astrometry of the MAD images of RX J1856. 5-3754 
was computed using as a reference the coordinates and 
position of 8 2MASS stars (Skrutskie et al. 2006), all 
non-saturated and evenly distributed in the field. Their 
pixel coordinates were measured by fitting their inten- 
sity profiles with a Gaussian function using the Graphical 
Astronomy and Image Analysis (GAIA) interfacqj. The 
coordinate transformation between the detector and the 
celestial reference frame was then computed using the 
Starlink package ASTROI^H. The rms of the astrometric so- 
lution was sa 0"17. After accounting for a 0"2 conservative 
astrometric accuracy of 2MASS, and the ~ 15 mas accu- 
racy (Skrutskie et al. 2006) on the tie to the International 
Celestial Reference Frame (ICRF), the overall error to be 
attached to our astrometry was finally 0"26. As a refer- 
ence for the RX J1856.5— 3754 coordinates at the epoch 

2 star-www.dur.ac.uk/ pdraper/gaia/gaia.html 
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of the MAD observation (2007.74) we used those of its op- 
tical counterpart reported in Kaplan et al. (2002), epoch 
1999.26, corrected for the measured proper motion. The 
uncertainty to be attached to the 2007.74 coordinates in- 
cludes the error on the proper motion extrapolation be- 
tween the two epochs, ~ 9 mas both in right ascension 
and declination, and the error on the reference coordi- 
nates which is 0'.'2 relative to the ICRF. We then com- 
puted the overall uncertainty on the position determina- 
tion of RX J1856. 5-3754 on the MAD images by adding 
in quadrature the error on our astrometric calibration and 
the uncertainty on its coordinates extrapolation, which is 
- 0"35. 

Since no 2MASS star is identified in the much less 
crowded RX J0420. 0—5022 field, we performed a two-step 
astrometric calibration using as a reference a grid of stars 
selected from the VLT H-band IS A A C image of the field 
(Lo Curto et al. 2007). By applying the same procedure 
as before, the astrometric calibration of the ISAAC im- 
age was performed using 8 2MASS stars as a reference, 
yielding an rms of ~ 0"22 on the astrometric solution. We 
then computed the coordinates of 7 field stars in common 
between the ISAAC and MAD fields, which made our sec- 
ondary astrometric grid. The astrometric solution yields 
an rms of ~ 0"09. By adding in quadrature all the errors 
of our astrometric chain, we then end up with an accu- 
racy of 0732 on our astrometry. As a reference for the RX 
J0420. 0—5022 coordinates, we used the Chandra coordi- 
nates of Haberl et al. (2004), which have a nominal error 
of 0"6 (90% confidence level). We accounted for the re- 
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cently measured proper motion upper limit (138 mas yr ; 
Motch et al. 2007) when we extrapolated the Chandra 
coordinates of RX J0420.0-5022 (epoch 2002.86) to the 
epoch of the MAD observation (2007.74). This yields a 
position uncertainty due to the proper motion of ~ 0'.'7. 
By adding in quadrature the error on our astrometric cal- 
ibration and the error on the original Chandra coordi- 
nates we ended up with an uncertainty of ~ 1" on the 
RX J0420. 0-5022 position at the epochs of the MAD ob- 
servations. 

The two panels in Fig. 1 show a zoom of the MAD im- 
age around the expected position of RX J 1856. 5— 3754 and 
RX J0420.0-5022. Despite of the sharper PSF of MAD, 
our images do not unveil any new object with respect to 
the ones identified in the long-exposure ISAAC H-band 
images of Lo Curto et al. (2007). In both cases, no can- 
didate counterpart is detected at the expected position, 
with the closest objects located several arcseconds away. 



3.2. Photometry 

To determine whether our observations can be used 
to better constrain the NIR emission properties of RX 
J1856. 5-3754 and RX J0420.0-5022, we estimated their 
K s -band brightness upper limits from the limiting mag- 
nitudes of the MAD images. To this aim, we have used 
as a reference the flux of a number of artificial stars sim- 
ulated with the template PSF measured on the images. 
Using different test values for the PSF flux normaliza- 
tion, artificial stars were simulated and added to the cor- 
responding images at different coordinates, randomly se- 
lected within a few arcsec from the target position. We 
then ran in loop the Daophot detection tool to deter- 
mine the PSF flux normalization corresponding to a 3<r 
detection limit. In all cases, aperture correction was ap- 
plied. We thus determined that 3cr detections would cor- 
respond to magnitudes of K s ~ 21.5 and K s ~ 20 for RX 
J0420.0-5022 and RX J1856. 5-3754, respectively. As an 
independent check, we have compared these values with 
the 3cr detection limits extrapolated from the fluxes of the 
faintest objects detected in our images. As a safe mea- 
sure, we used the ISAAC H-band images of Lo Curto et 
al. (2007) to confirm that they were indeed genuine de- 
tections and not background fluctuations. Following this 
empirical approach, we derived detection limits consistent 
with the previous ones. The different limiting magnitudes 
for the two images, despite comparable integration times, 
has to be ascribed mainly to the different image qual- 
ity between the two frames, being much better for RX 
J0420.0-5022 than for RX J1856.5-3754 (tf'28 with re- 
spect to 0'.'40, see Sect. 2.2). This is indeed a significant 
factor when the signal-to-noise ratio is background lim- 
ited. Moreover, the RX J1856.5— 3754 image is affected by 
a higher background level which is ~ 25 % higher around 
the target position. 



4. Discussion 

Our measured K s -band spectral flux upper limits are well 
above the extrapolation of the optical spectrum of RX 
J1856.5-3754 (van Kerkwijk & Kulkarni 2001), which 
would predict K s ~ 26.5, and the Ray leigh- Jeans extrap- 
olation in the optical/NIR domain of the XMM-Newton 
spectrum of RX J0420.0-5022 (Haberl et al. 2004). We 
note that, for RX J1856.5— 3754, the K s -band spectral flux 
upper limit is shallower than the H-band one of Lo Curto 
et al. (2007), while for RX J0420.0-5022 it is slightly 
deeper. Thus, we are not yet able to constrain the NIR 
spectrum of RX J1856. 5-3754 and RX J0420.0-5022. 

A blackbody spectrum produced by the neutron star 
surface is obviously a possibility. It is clear that for the RX 
J1856.5-3754 and RX J0420.0-5022 distances of - 160 
pc (van Kerkwijk & Kaplan 2007) and ~ 350 pc (Posselt 
et al. 2007), respectively, and for any reasonable combi- 
nation of temperature and emitting area, the predicted 
flux would fall well below our upper limits. A blackbody 
spectrum produced by a fallback disk (Perna et al. 2000) 
could yield to a flux compatible with the present NIR flux 
upper limits of RX J1856. 5-3754 and RX J0420.0-5022, 
depending on the actual disk size and accretion rate (see 
also Lo Curto et al. 2007). A power-law spectrum could 
be produced by non-thermal emission from the neutron 
star magnetosphere. In FigJ5] (left), we plotted the NIR 
luminosity (L^ir) vs the rotational energy loss (E) for 
different classes of INSs: rotation-powered neutron stars, 
magnetars, and XDINSs. Our upper limits do not rule 
out a NIR emission efficiency, t]nir = Lnir/E, compara- 
ble to that of rotation-powered neutron stars, or higher. 
This would be possible, for instance, if the fraction of en- 
ergy emitted in the relativistic wind is much less than in 
rotation-powered neutron stars or if the NIR emission is 
powered not by the neutron star rotation but, e.g., by its 
magnetic field. Although XDINSs are less luminous than 
magnetars in absolute terms (Fig[2J left), they might have 
a comparable NIR emission efficiency. This is shown in 
FigE] (right) where the NIR emission efficiency is plot- 
ted vs. the dipole magnetic field B inferred from the spin 
down. We warn here that NIR observations of the radio 
pulsar PSR J1119-6127 (Mignani et al. 2007b) suggest 
a NIR emission efficiency comparable to that of other 
rotation-powered neutron stars, despite of its relatively 
high magnetic field (B ~ 4.1 x 10 13 G). If the inferred 
magnetic field value of PSR J1119— 6127 is correct, this 
might argue against a connection between the NIR emis- 
sion efficiency and the magnetic field. 

Much deeper observations are required to better con- 
strain the XDINS optical/NIR emission properties and 
to investigate their possible connection with the magne- 
tars. In this respect, particularly important would be the 
study of 1RXS J214303.7+065419, possibly the XDINS 
with the highest magnetic field, as hinted by the detec- 
tion of the X-ray absorption feature (Zane et al. 2005). Its 
anomalously high optical emission is indeed incompatible 
with both thermal emission from the neutron star surface 
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Fig. 2. (left) INS NIR luminosities vs. rotational energy loss E. Filled circles, squares and triangles correspond to 
rotation-powered neutron stars, magnetars, and XDINSs, respectively. Plots are updated from Mignani et al. (2007b) 
and Lo Curto et al. (2007), to which we refer for the complete source list and for the timing parameters. Luminosities 
are computed in the K s band for magnetars and for rotation-powered neutron stars. For the XDINSs, we have plotted 
the luminosity upper limits computed from the H-band observations of Lo Curto et al. (2007). For RX J0420. 0—5022 
only, we have plotted the slightly deeper luminosity upper limit obtained from our K s -band observations. The red 
dashed line is the linear fit to rotation-powered neutron stars, while the solid line marks the limit case Log(LNjn) = 
Log(-E'). (right) INS NIR efficiencies, "q^iR = Lnir/E, vs. dipole magnetic field B inferred from the spin down. Note 
that for those XDINSs for which only an upper limit on E is available, including RX J0420. 0—5022, the NIR emission 
efficiency is unconstrained, as shown by the vertical dashed lines. 



and rotation-powered non-thermal emission (Zane et al. 
2008). 
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